Fluctuation-Induced Casimir Forces in Granular Fluids

We have numerically investigated the behavior of driven non-cohesive granular media and found that two fixed large intruder particles, immersed in a sea of small particles, experience, in addition to a short range depletion force, a long range repulsive force. The observed long range interaction is fluctuation-induced and we propose a mechanism similar to the Casimir effect that generates it: the hydrodynamic fluctuations are geometrically confined between the intruders, producing an unbalanced renormalized pressure. An estimation based on computing the possible Fourier modes explains the repulsive force and is in qualitative agreement with the simulations.Comment: 4 pages, 3 figures. Accepted in Phys. Rev. Letter

A Particle-based Multiscale Solver for Compressible Liquid-Vapor Flow

To describe complex flow systems accurately, it is in many cases important to account for the properties of fluid flows on a microscopic scale. In this work, we focus on the description of liquid-vapor flow with a sharp interface between the phases. The local phase dynamics at the interface can be interpreted as a Riemann problem for which we develop a multiscale solver in the spirit of the heterogeneous multiscale method, using a particle-based microscale model to augment the macroscopic two-phase flow system. The application of a microscale model makes it possible to use the intrinsic properties of the fluid at the microscale, instead of formulating (ad-hoc) constitutive relations

Optimized Verlet-like algorithms for molecular dynamics simulations

New explicit velocity- and position-Verlet-like algorithms of the second order are proposed to integrate the equations of motion in many-body systems. The algorithms are derived on the basis of an extended decomposition scheme at the presence of a free parameter. The nonzero value for this parameter is obtained by reducing the influence of truncated terms to a minimum. As a result, the new algorithms appear to be more efficient than the original Verlet versions which correspond to a particular case when the introduced parameter is equal to zero. Like the original versions, the proposed counterparts are symplectic and time reversible, but lead to an improved accuracy in the generated solutions at the same overall computational costs. The advantages of the new algorithms are demonstrated in molecular dynamics simulations of a Lennard-Jones fluid.Comment: 5 pages, 2 figures; submitted to Phys. Rev.

Temperature scaling in a dense vibro-fluidised granular material

The leading order "temperature" of a dense two dimensional granular material fluidised by external vibrations is determined. An asymptotic solution is obtained where the particles are considered to be elastic in the leading approximation. The velocity distribution is a Maxwell-Boltzmann distribution in the leading approximation. The density profile is determined by solving the momentum balance equation in the vertical direction, where the relation between the pressure and density is provided by the virial equation of state. The predictions of the present analysis show good agreement with simulation results at higher densities where theories for a dilute vibrated granular material, with the pressure-density relation provided by the ideal gas law, are in error. The theory also predicts the scaling relations of the total dissipation in the bed reported by McNamara and Luding (PRE v 58, p 813).Comment: ReVTeX (psfrag), 5 pages, 5 figures, Submitted to PR

A Combined Spectrophotometer and Fluorometer to Demonstrate the Principles of Absorption Spectroscopy

A dual-function student-crafted instrument is described as part of a laboratory activity aimed to teach both the principles and practical aspects of absorption spectroscopy to secondary and introductory undergraduate students. Using minimal changes in an arrangement that is based on interlocking bricks and low-cost components, both a fluorometer and photometer have been constructed. The former demonstrates the principles of the Beer–Lambert law visually and quantitatively by acquiring the spatial light attenuation through a fluorescent sample. The latter then demonstrates its practical application in a visible-light spectrometer by measuring the absorption spectrum of an aqueous permanganate solution

Algorithm for numerical integration of the rigid-body equations of motion

A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in terms of either principal axes or quaternions. Due to specific features of the algorithm, orthonormality and unit norms of the orientational variables are integrals of motion, despite an approximate character of the produced trajectories. It is shown that the method presented appears to be the most efficient among all known algorithms of such a kind.Comment: 4 pages, 1 figur

Quasi-classical Molecular Dynamics Simulations of the Electron Gas: Dynamic properties

Results of quasi-classical molecular dynamics simulations of the quantum electron gas are reported. Quantum effects corresponding to the Pauli and the Heisenberg principle are modeled by an effective momentum-dependent Hamiltonian. The velocity autocorrelation functions and the dynamic structure factors have been computed. A comparison with theoretical predictions was performed.Comment: 8 figure

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.Comment: EASC 2014 conference proceedin

Diffusive Spreading of Chainlike Molecules on Surfaces

We study the diffusion and submonolayer spreading of chainlike molecules on surfaces. Using the fluctuating bond model we extract the collective and tracer diffusion coefficients D_c and D_t with a variety of methods. We show that D_c(theta) has unusual behavior as a function of the coverage theta. It first increases but after a maximum goes to zero as theta go to one. We show that the increase is due to entropic repulsion that leads to steep density profiles for spreading droplets seen in experiments. We also develop an analytic model for D_c(theta) which agrees well with the simulations.Comment: 3 pages, RevTeX, 4 postscript figures, to appear in Phys. Rev. Letters (1996

Thermal conductivity of the Toda lattice with conservative noise

We study the thermal conductivity of the one dimensional Toda lattice perturbed by a stochastic dynamics preserving energy and momentum. The strength of the stochastic noise is controlled by a parameter $\gamma$. We show that heat transport is anomalous, and that the thermal conductivity diverges with the length $n$ of the chain according to $\kappa(n) \sim n^\alpha$, with $0 < \alpha \leq 1/2$. In particular, the ballistic heat conduction of the unperturbed Toda chain is destroyed. Besides, the exponent $\alpha$ of the divergence depends on $\gamma$